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Abstract 

An experimental technique has been developed for the measurement 
of post-shock temperatures in a wide variety of materials, including 
those of geophysical interest such as silicates. The technique uses 
•an infra-red radiation detector to determine the brightness temperature 
of samples shocked to pressures in the range 5 to ^30 GPa; in these 
experiments measurements have been made in two wavelength ranges 
(4.5 to 5.75p; 7 to 14p). Reproducible results, with the temperatures 
in the two wavelength bands generally in excellent agreement, have 
been obtained for aluminium-2024 (10.5 to 33 GPa; 125 to 260°C), stain- 
less steel-304 (11.5 to 50 GPa; 80 to 35Q°C), crystalline quartz 
(5.0 to 21.5 GPa; 80 to 250°C), forsterite (7.5 to 28.0 GPa; v30 to 
160°C) and Bamble bronzite (6.0 to 26.0 GPa; ^30 to 225°G). 

These results are generally much higher at low pressures (where 
they may even be in excess of the calculated shock temperatures) than 
the values calculated assuming a hydrodynamic rheology and isera tropic 
release parallel to the Hugoniot but tend towards them at higher 
pressures. In aluminium- 20 24, the theoretical post-shock temperatures, 
assuming a fluid- like rheology, are 35 to 218°C, for the pressure range 
10.5 to 33 GPa. However, the results are in considerably better 
agreement with values calculated assuming elasto-plastic behaviour 
(80 to 270 °C) which probably also causes the high measured temperatures 
for stainless steel. In forsterite the measured values ranged from 
65°C to 9.6 GPa (there was no detectable rise at 7.5 GPa) to 156° at 
28.0 GPa, whereas the "hydrodynamic values" were 30 to 120°C. Values 
obtained for quartz were in excellent agreement with those calculated 
by Mashlmo et al. (1979) using release adiabat data. It is concluded 
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that release adiabat data should be used, wherever available, for 
calculations of residual temperature, and that adequate descriptions of 
the shock and release processes in minerals need to be more complex than 
generally assumed. 

Introduction 

Shock wave equation of state data have long been used in the 
interpretation of impact metamorphism (e.g., Stbffler, 1971, 1972) and 
of density depth profiles of the earth obtained from seismic data 
(e.g., Al'tschuler, 1965). However, one of the limitations to these 
uses of Hugoniot data is the uncertainty in the temperatures reached 
both during the passage of the shock wave through the material and 
after unloading. 

Shock temperatures have previously been calculated using Hugoniot 
data and the Mie-Griineisen theory (e.g., Walsh and Christian, 1955; 

Wackerle, 1962; Ahrens et al. , 1969, McQueen et al . , 1970) which should 
yield reliable results for metals provided the rheology corresponds to 
fluid- like behaviour. However, the Mie-Griineisen theory, when the 
material is treated as a single phase, is Inadequate over much of the 
pressure range for which shock wave data exist for silicates, since these 
all undergo major phase changes which may involve substantial changes in 
thermodynamic properties and internal energy. Calculated Hugoniot temperatures 
can then be used as a starting point to obtain post-shock temperatures assuming 
isentropic release. Unfortunately, there are few cases where there are 
sufficient data to either experimentally determine or theoretically 
calculate isentropic release paths; moreover the assumption of isentropic 
unloading from shock states has never been explicitly verified 
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experimentally. Consequently, even if the Hngoniot temperatures were 
calculated correctly, large uncertainties could still exist in residual 
temperatures . 

Calculations of post-shock temperatures in silicates, assuming 
release along isentropes lying above the Hugoniot (Fig. 1), with the 
release volume greater than the initial volume, lead to values that 
appear too low to account for some of the effects seen in recovery 
experiments (eg. Gibbons, 1974). "Where release adiabat data exist for 
silicates, they have been used in the calculation of post-shock 
temperatures (eg. Gibbons and Ahrens, 1971; Ahrens and O’Keefe, 1972). 
This always leads to much higher, and possibly more credible, release 
temperatures (as is shown by the comparison in Table 1) largely because 
the release paths lie below the Hugoniot. 

Early experiments to measure post-shock temperatures were limited 
to metals; using a photo-multiplier tube Taylor (1963) found values for 
copper shocked to pressures in the range 90 to 170 GFa that were in good 
agreement with those calculated by McQueen and Marsh (1970) using Mie- 
Grilneisen theory and fluid rheology, although later work by King et al. , 
(1967) and Von Holle and Trimble (1976) using infra-red radiation 
detectors revealed considerable discrepancies between observed and 
calculated temperatures for copper shocked to pressures below 80 GPa. 
However, no attempts were made to extend this type of experiment to 
non-metals, largely because the supposedly lower temperatures, coupled 
with the low sensitivity of available detectors, meant that experiments 
such as those of King et al. wore not feasible. Recent improvements in 
detector technology have now made It possible to design a system capable 
of measuring post-shock temperatures in silicates; the availability of 
such data should help resolve the current uncertainties. 
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Experimental Technique 

Initial experiments were designed with the aim of investigating 
residual temperatures in silicates of geophysical interest shacked to 
pressures up to 30 GPa. The materials chosen were natural crystal 
quarts (cut perpendicular to the c axis) Bamhle bronzite (from Bamble, 

Norway) , which consists of large natural single crystals closely described 
by the formula (Mg^ g^-Fe^ 14^^3 k ave a porosity of 1 to 3 %, and 

synthetic crystal forsterite (cut perpendicular to the c-axis). For 
completeness, the materials used as driver plates in the experiments, 
namely aluminium- 2024 and stainless steel-304, were also studied. 

The method developed involves the monitoring of radiation from the back (free) 
surface of a shocked sample with an infra-red radiation detector whose 
output may rh en be used to determine the brightness temperature of the 
sample. Since the residual temperatures for the pressure range to be 
investigated were expected to be of the order of 400°K, and the Planck 
distribution law gives a maximum in spectral radiance between 3 and 9 V 
for black bodies radiating at temperatures from 1000 to 300°K (Touloukian 
and DeWitt, 1972) infra-red detectors were a logical choice for this study. 

In addition, since silicates behave as fairly good black bodies in the 
infra-red beyond ^5y, with the exception of a minimum in emissivity at 
^9y (Lyon, 1965), the radiative output of the sample is also maximised. 

A schematic plan of the experimental configuration is shown in Fig. 2. 

a) Production of the shocked state. 

In these experiments, the shocked state was produced in the sample 
by the impact of a gun launched flyer plate. The technique is described 
in detail by Gibbons (1974) and Ralkes (1978) and has been used by a 
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number of authors (eg. Ahrens et al., 1971; Ahrens and Gaffney, 1971; 

King and Ahrens, 1976). The velocity of the projectile is determined 
from the interval between obscurations (by the front of the projectile) 
of the two laser beams (Fig. 2); typical uncertainties in velocity are 
±0.05 km/s. The shock pressure in the sample is then calculated using 
the impedance match method (eg. Duvall and Fowles, 1963) and the known 
Hugoniots of flyer plate, driver plate and sample materials. Estimates 
of the uncertainty in pressure determination, based on the accuracy of 
velocity measurement and scatter of the Hugoniot points, are ±0.5 to 
1.0 GPa in the silicates. 

After passing through the mylar window sealing the end of the 
barrel, the projectile makes contact across a shorting target made from 
thin copper foil and having a negligible effect on the projectile velocity, 
and finally impacts the target assembly, (Fig. 2) . In order to achieve 
nearly one- dimensional planar flow upon impact of the projectile with 
the target, the target is carefully aligned using the normal reflection 
from the target surface of a laser beam shining down the centre of the 
barrel. The target assembly consists of a driver plate (1.5mm aluminium- 
2024 or stainless steel- 304), followed by a 3mm thick silicate sample, 

15mm in diameter; in the case of shots on metals the driver plate, now 
3mm thick, is the sample. 

b. Temperature measurement 

The brightness temperature of the back (free) surface of the sample 
is determined from the output from au infra-red detector. This detector 
is mounted above the impact chamber (Figure 2) , and monitors the back 
face of the sample via a mirror and optical system. The latter ensures 
that only the centre (^0.8 cm dia) of the sample is viewed, reducing the 
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contribution from radial release waves (also somewhat lessened by 
using a circular sample) and increases the efficiency of the detection 
system. The detector is connected via an amplifier to two oscilloscopes. 
One is triggered by the passage of the projectile past the first laser 
beam of the timing system, and records the detector output at a rate 
of 50 ys/div. This provides a back-up record in case of failure of the 
higher time-resolution recording, and a means of checking that no 
temperature signals are generated prior to the passage of the shock wave 
through the sample assembly. (Note that it also provides another means 
of determining the projectile velocity. ) The second ocilloscope is 
triggered by the contact of the flyer plate with the shorting target, 
which is approximately 15 mm in front of the driver plate, just prior 
to impact; this writes at 5 ys/div, and it is the primary record that 
is used in temperature determination. The interpretation of this 
record is based on the assumption that the oscilloscope sweep is triggered 
at the moment of contact of the flyer plate with the shorting target; 
this was checked by analysing the timing of signals seen on the back-up 
record (50 ys/div) and by using a delayed trigger on the oscilloscope, 
and appears valid to within ±0.5 ys, which is close to the rise time of 
the systems. A typical record (such as that for hronaite at 15.5 GPa, 

Pig. 8c) shows a sharp rise in signal level as the shock wave arrives 
at the free surface of the sample followed by a level portion corresponding 
to the residual temperature, and then a subsequent rise due to air shocks 
generated at the end of the sample chamber and the destruction of the 
mirror. At high pressures, especially for metals which have a second 
peak in the record (Fig. 7), the level portion may not be clearly visible; 
in this case the post-shock temperature is determined from the signal 
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0.5 (is (InSb) or 0.75 ps (HgCdTe) after the free surface arrival. (These 
values are the rise times of the detection systems.) Problems in resolving 
the residual temperature were also encountered for quartz where an initial 
flash (Fig. 8a) made identification of the level portion uncertain at the 
highest pressures. 

The InSb and HgCdTe detectors were operated at 77 °K. 

They are enclosed in dewars cooled by liquid nitrogen; 
both were obtained from the Santa Barbara Research Centre (Goleta, 
California) and typical response curves, obtained from S.B.R.C., are 
shown in Fig. 3. Additional operational details are listed in Appendix A. 
A filter was used to limit the bandwidth of the InSb detector to 4.5 to 
5.75p in order to minimise the possibility of radiation from the metal 
driver plate, or metal-sample interface, being transmitted through the 
sample and causing errors in the temperature determination. (As can 
be seen from the transmission scans of Fig. 4, which were obtained using 
a Perkin- Elmer Model 180 Infra-red Spectro-photometer in a manner 
analogous to that described in Burns (1970) , this is unlikely to be a 
problem for quartz or bronzite but the forsterite does transmit 
significantly up to 5p.) 

The InSb detector was used with a variable gain amplifier (Raikes, 
1978); the system rise time, which is essentially limited to v0.2 ps 
by the chip itself, was ^0.5 ps depending on the gain setting used. 
Although the HgCdTe is a faster material (<0.1ps), the rise time of the 
detector-amplifier system was ^,0.8ps because an electronic filter with 
a high frequency cut-off of 5MHz had to be used to reduce the large 
amount of high frequency noise which would otherwise have made accurate 
measurement of temperature impossible. 
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In order to convert the voltage record into a brightness temperature, 
the detector must be calibrated. This is best performed by heating the 
sample -£fi &1&L to a known temperature, and recording the corresponding 
voltage output of the detector-amplifier system. This is easily done 
for metals, but would be very hard for the non-metals studied since they 
are extremely brittle and hard to heat in the experimental configuration 
without cracking. Instead, initial estimates of temperature were based 
on the assumption that the silicates behaved as black bodies using a 
calibration curve obtained for a "black” body, graphite. For quartz in 
range 5 to 8p this is a reasonable approximation since the emissivity is 
greater than 0.9 (Touloukian and DeWitt, 1972), and is probably justified 
for the Bamble bronzite in the InSb range, since it does not transmit. 
However, the forsterite has a 20% transmittance at 4.5p (dropping 
rapidly to less than 5% at Sp), and all the silicates studied have a 
large drop in emissivity at around 9p. This is clearly visible in the 
comparison of black body and quartz emittance spectra at temperatures 
from 250 to 500°K presented by Lyon (1965). Corrections for sample 
emissivity were determined using available data, which are unfortunately, 
largely for powdered minerals or rocks. The values of the emissivity for 
quartz and dunite (primarily forsterite) given by Buettner and Kern (1965) , 
Lyon (1965) and Touloukian and DeWitt (1972) were used to calculate 
correction factors as follows: 

The detector output S for an operating wavelength range from Yx 
Y 2 P may be expressed as 



E(X,T)D(X)P(X,T)dX 


(D 
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where E (A,T) - E(A) - emissivity (assumed independent of T) 

D(X) = detector response, P(A,T) = Planck's function 

A = wavelength, X = absolute temperature 
This may be integrated numerically and used to derive the ratio of the 
signal obtained for a silicate at temperature T to that for a black body 
at the same temperature. The correction factors are listed in Table 2 
the values for bronzite were estimated from those for quartz and dunite 
based on a comparison of absorption spectra of the three materials. 

Since the post-shock temperatures are measured after the interaction 
of the shock-wave with the free surface, and that interaction will cause 
roughening of the surface, a correction factor derived for a rough surface 
might be more appropriate; however, the values listed in Table 2 for 
polished surfaces will be used since these should yield an upper bound 
on the temperature. Two additional factors should be taken into account: 
one is that the absorption peak may shift during shock compression 
(e.g., Goto et al. , 1979), and broadening of the absorption bands for 
Si 02 has "been observed in samples recovered after shock compression to 
pressures up to 52 GPa (Mashimo et al . , 1979; Stoffler, 1974). The 
second is the possibility of tribo luminescence, or some other form of 
non- equilibrium radiation such as might be associated with a phase 
change; in these cases the emlssivity may even exceed unity. Because 
of these uncertainties, the black body temperatures may well be more 
reasonable estimates of the residual values and will be used later in a 
comparison of observed and calculated temperatures; the corrected values 
probably represent upper bounds to the post-shock temperatures. 

Typical calibration curves are plotted in Pigs. 5 and 6, and 
Table 3 gives the power law fits to the curves. Both detectors gave 
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extremely reproducible calibration curves , as is demonstrated by the 
two sets of points for aluminium in the case of XnSb (these were obtained 
several weeks apart with a number of shots fired in between) . 

Sources of error. 

Errors in temperature measurement can basically arise from two 
causes — • those related to the sample, and those originating from out- 
side sources. Prime among the latter is contamination of the signal by 
radiation from air shocks which can he of extremely high temperature. 

Owing to the geometry of the apparatus these should not be important 
prior to the arrival of the shock wave at the free surface of the 
sample, and, indeed, no earlier signal rises were detected. The only 
air shock likely to affect the post-shock temperature measurement would 
be one generated at the back of the sample itself, which should be 
eliminated by the sample vacuum chamber that is pumped down to ^5m Torr 
(£0.7 Pa). To reduce further possible radiation from residual gases 
within this chamber heated by compression due to the shock wave, the 
chamber was flushed oat with helium prior to each shot. Radiation 
from later air shocks, such as that generated at the end of the sample 
chamber as the window breaks due to the shock wave travelling down the 
chamber walls and impact of ejecta from the samples, is clearly visible 
on each record, and ultimately causes the detection system to saturate. 

If the samples were transparent, then radiation from the metal driver 
plate could add to the signal, but the rise should then precede the 
free surface arrival of the shock wave, and this does not appear to be 
the case. Anyway, the emissivity of the metal is substantially lower than 
that of the sample, so this effect should be relatively small, and both 
detectors were chosen to operate at wavelengths where the samples are 


nearly opaque. (Forsterite does have a 2Q% transmittance at 4.5p 
dropping rapidly to less than 5 % at 5p, and so may show minor effects 
due to transmitted radiation for the InSb detector.) Both the driver 
plate and the sample surface in contact with it were polished prior 
to mounting of the sample in order to minimise the "porous" surface 
interactions that could give rise to considerable heating (see e.g,, 

Urfciew and Grover, 1974); this also reduces the likelihood of air 
being trapped in this interface, but this should be removed by the 
evacuation of the sample chamber. In order for the signal to be 
contaminated by radiation from the metal or the driver-sample interface, 
large changes must take place in the transmissivity of the sample under 
shock conditions; although changes have been reported in sapphire 
(Urtiew, 1974), these were decreases and at much higher pressures. 

It seems that this is not a likely source of error. Vibration of the 
detector- amplifier system could conceivably affect the output. However, 
the detection system was physically clamped in isolation from the gun, 
and vibrational effects were not apparent except for some HgCdTe shots 
where a negative signal of short duration (vLO us) , obviously non- thermal 
in origin, was observed prior to impact with the shorting target, but 
tbe output returned to the aero level before the free surface arrival. 

Two main sources of errors associated with the sample behaviour under 
shock are changes in emissivity and the effects of non-uniform heating. 

The former may be investigated by comparing the brightness temperatures 
obtained at different wavelengths, since the change would probably not 
be constant as a function of wavelength. Changes in emissivity may be 
related to phase changes, changes in surface properties and triboluminescence. 
In a recent work on residual temperatures in copper. Von Holle and Trimble 
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(1976) determined temperatures "by using the ratio of detector signals 
obtained at different wavelengths. They believed this would reduce 
the likelihood of error due to changes in emissivity and the effects 
of surface processes, and yield a relatively unbiassed estimate of the 
post-shock temperature. Unfortunately, the signal to noise ratio is 
low, as at low temperatures and especially for metals where the emissivity 
is low, small errors in the measurement of detector output can lead to 
large changes in the ratio of the signals from the two detectors and 
totally unreasonable temperatures. In these experiments the ratio 
technique did not prove useful, although for the higher pressure shots 
the calculated ratio temperatures were, in fact, in good agreement with 
the separate brightness temperatures. 

The effects of non-uniform heating may he more severe; they were 
estimated by calculating the temperature that would be measured if 10% 
of the surface were 100 or 2Q0°K hotter than the remaining 90%. For 
a 100° excess, the discrepancy between "measured" and mean temperatures 
is close to the accuracy of measurement, i.e ^±10%. However, for a 
200° excess the differences may reach v30° (with the measured temperature 
being an over-estimate), although such a temperature distribution seems 
rather unlikely because of the energy partition that would be required. 
Non-uniform shock heating has been clearly demonstrated in silicates 
where such features as adiabatic shear zones (e.g., Grady, 1977) 
deformational twin lamellae and localised production of glass (e.g., 
in bronzite. Gibbons, 1974) have been observed. This non-uniform heating 
is believed to contribute to the complete loss of strength once the 
Hugoniot elastic limit is exceeded. The effect of surface processes 
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such as jetting can also bias the temperature measurements; these will 

presumably be more important in metals where the optical depth is of the 

order of angstroms, than in silicates where it is microns. To reduce the 

likelihood of jetting, the sample surfaces were polished, but not to a 

high gloss as this would reduce the emissivity. Surface processes are 

probably more important in metals where the infra-red optical depths 
-10 ' 

are MX) m, than silicates where they are microns. 

Observations 

a. Stainless steel-304 and Aluminium- 2024 

Typical detector output records are shown in Tig. 7: these are for 
the InSb detector, but records obtained with the HgCdTe detector were 
essentially similar except for the longer rise time and the fact that 
the amplifier saturated prior to the peak T 2 » The origin of this later 
peak is unclear: the temperature correlates well with pressure, shock 
and free surface velocities for both stainless steel and aluminium, 
the correlation curves for the two materials are distinct even if a 
common temperature calibration curve is used it appears to be a material 
property, rather than some effect common to all shots, such as the 
compression of residual gas within the sample chamber (it was not observed 
in shots on silicates at similar pressures) and may be related to 
localised heating on break up of the sample or to the action of radial 
release waves. 

The temperatures determined are listed in Tables 4 and 5: agreement 
between the two wavelength ranges is good, with the difference not ex- 
ceeding the estimated uncertainty in measurement. Also tabulated are 
the corresponding black body temperatures, which should place a lower 
bound on the residual temperatures; the values for aluminium are probably 
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too low, since the maximum reported emissivity for roughened oxidised 
aluminium- 20 24 is ^.4 (Touloukian and Ho, 1972) which would imply 
temperatures of 90° to 26Q°C. 
b. Silicates 

Typical records obtained using the InSb and HgGdTe detectors are 
shown in Figs. 8 and 9. Unlike quartz and bronzite, forsterite has 
very different records for the two wavelength ranges, probably because 
it transmits up to 20% in the InSb range; the InSb temperatures may 
thus be contaminated by radiation from the sample- driver interface. 

The records for quartz have a characteristic short duration 'flash 1 
which is probably associated with tribo luminescence, a phenomenon 
documented in quartz shocked to similar pressures by Nielson et al. , 
(1961); the black-body temperatures corresponding tf‘< this peak are 
unlikely to be significant. Black body and corrected temperatures 
for the silicates studied are listed in Tables 6 to 8. 


Comparison of Observed and Calculated Temperatures 
a. Calculational techniques 

The most widely used method of calculating shock, and hence post- 
shock, temperatures is probably that developed for metals by Walsh and 
Christian (1955) who derived a differential equation relating the shock 
temperature Tjj to the volume on the Hugoniot and material properties. 
This equation has the general solution 


\ ( V - *« •* < - f 


b(v) dv) + 


( 2 ) 


V ,H 

exp ( - f b(v) dv) f f (v) exp ( f b(v)dv) dv 

V 0 V v ^0 


Here T and v are the initial temperature and volume, C v the specific 



heat at constant volume, b = \fiy = v w * iei:e Y “ Gruneisen's 
parameter and f (v) ~ . ( v o~ v ^ ~ evaluated at pressure, volume, points 


(P,v), along the Hugoniot. 

A common assumption is that y has the form 

Y " T °tF 

where n £ 1 if n = 1, then b is constant and (2) reduces to the form 
actually derived by Walsh and Christian. Equation (2) is reduced to 
a difference equation and solved iteratively along the Hugoniot for 

%• 

The derivation of (2) is based on two main assumptions: firstly 

that an ordinary fluid type equation of state is valid, which ignores 

the effects of rigidity or alas to-p las tic work, and secondly that 

thermodynamic equilibrium exists in states behind the shock front. 

In addition, it is strictly valid only where the shock state is reached 

by a single step and not in the two wave- region associated with the 

Hugoniot elastic limit or phase changes; it must thus be modified for 

application to silicates such as those studied here. In the case of 

quartz, Wackerle (1962) modified the approach of Walsh and Christian by 

calculating temperatures along an "equilibrium" Hugoniot obtained from 

* 

a segmented linear fit to the plot of effective shock velocity U against 
effective particle velocity u*, where U and u* are given by 


u * = W <v o - V 1 1/2 


(3) 


“* - IV-o - V ]1/2 

ire v„ is the volume corresponding to a shock 
ti 

iese represent the true velocities only at high enough 


pressure P^. 


( 4 ) 


the two-wave structure is overdriven. ) Such equ i librium Hugoniots may 
be defined for quartz anf forsterlte using Wackerle’s 1962 data for 
quartz and the single crystal datum reported in Appendix B, the data 
for poly crystalline forsterite (McQueen, 1968) and ^490 porous poly- 
crystalline forsterite (Ahrens et al. , 1971). However, the data for 
bronzite are so scattered probably because of varying porosity, that 
no satisfactory equilibrium Hugoniot can be derived; instead, Walsh 
and Christian temperatures were calculated along Hugoniots derived from 

Birch Murnaghan adiabats using tlie constants of Table 9 and an initial 
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specific volume of 0.304 cm. /g (for details of this method see, e.g., 
Davies, 1974). A theoretical Hugoniot was also used in calculations 
involving forsterite. 

An alternative approach was described in detail by Ahrens efc al. 
(1969) . The increase in internal energy AE^ of a material shocked to 
a state with volume v^ and pressure P^ is given by the Rankine- Hugoniot 
conservation equation increase in internal energy resulting from iso- 
thermal compression at T q from an initial volume to a final volume V H 
plus isovolumie heating to the shock temperature T^. In this case 

AE h - 1/2 CE H + P 0 ) CV 0 - V H ) = <?$*- - P) T dv + 

v o 

where P,^ is the pressure on the isotherm, at volume, V^. 

This is solved for P„, and T„ determined from the expression for 

1 Ji 

the pressure difference between the Hugoniot and the isotherm 

P H = P T - 1. ? C v dT 

T 

o 

where b = y/v is assumed constant. 
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This formulation is easily adapted for inclusion of the Hugoniot 
elastic limit, and for calculating shock temperatures in the high 
pressure regime for materials -which undergo phase changes. Both quartz 
and hronzite do so in the pressure range investigated (starting at 
15 GPa and ^ 13.5 GPa respectively) hut neither is complete until 
£ 40 GPa. The highest pressure measurements here are thus in the 
mixed phase regime - in quartz at 20 GPa the phase change is $ 25 % 
complete (Grady et al. , 1974). Calculations in the mixed phase regime 
are extremely complicated, the Hugoniot state may be assumed to be 
a mixture of high and low pressure phases in thermal and mechanical 
equilibrium and equations derived which can be solved numerically for 
the mass fraction of transformed material and the Hugoniot temperature 
at a series of points on the mixed phase Hugoniot (Ahrens et al. , 1969) - 
However, in view of the uncertainties in the mechanism of the phase 
change, and in bronzite its nature (the high pressure phase is probably 
majorite (Smith and Mason, 1970), but may be ringwoodite plus stishovite 
(Ahrens and Gaffney, 1971)) and the pressure range studied, the effects 
of the phase change were ignored, except in so far as they changed the 
shape of the equilibrium Hugoniot for quartz. The calculation of post- 
shock temperatures in the mixed phase regime is further complicated by 
the behaviour of the high pressure phase an release, when it may transform 
to an amorphous form (Grady, et al. , 1974). 

The calculational methods derived by Walsh and Christian and 
Ahrens et al. both yield values for Hugoniot temperatures, and the residual 
temperatures must he determined from these. It is generally assumed that 
tbe release is isentropic, in which case the residual temperature Is 


18 


given by 



(7) 


where V_ — volume on release, often assumed equal to V . For eases where 
it o 

the release path is known, the residual temperature may he calculated 
directly, as described by Gibbons (1974) : the energy in the Hugoniot state 
is equated to the change in internal energy due to the rise in temperature 
from the initial value T to the residual value TL plus the energy 

O K * 

change on release. Thus 




1/2 CP H + Po> P 0 - V - / ^release I C p dl 
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( 8 ) 


H 




where V B is the volume on release and C the specific heat at constant 
K P 

pressure. 


b. Results for Aluminium and Stainless Steel 

Since the Walsh and Christian approach was developed specifically 
for application to metals (in the absence of phase changes) where the 
Hugoniot elastic limits are low (£l.O GPa) , the application of this 
technique should yield results in good agreement with the experimental 
observations. Fig. 10 (a) and (b) show the values of shock (dashed 
lines) and post-shock temperature (solid curves) calculated for stainless 
steel-304 and aluminium- 20 24 using this approach; they are the same as 
those given by McQueen et al. , (1970) . Also plotted are the observed 
values, and, as can be seen, there is practically no agreement. In 
fact, the measured residual temperatures are, at low pressures, in excess 
of the calculated Hugoniot temperatures. Discrepancies between observed 
and theoretical residual temperatures in metals have also been reported 
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by other workers (e,g., Von Holle and Trimble, 1976) using similar 
experimental techniques, and the question of the validity of the 
measurement arises. Certainly, the measured values may be too high 
because of surface processes or changes in emissivity — though they 
are consistent in two wavelength ranges — but the fact that they 
tend towards the theoretical values with increasing pressure suggests 
that the theory does not include some effect that dominates at low 
pressures. (Allowing for reasonable emissivity changes, the observed 
values are, in fact, still higher than those calculated.) The observed 
temperatures for steel may at first seem high when compared to the 
small amount of heating apparently observed in steel containers used in 
recovery experiments; however, these containers are not examined 
immediately, and the initial post-shock temperatures will quickly decay 
owing to thermal conduction. This is borne out by the observations of 
Schneider and Stilp (1977) who used thermocouples to measure the 
temperatures within large steel targets as a function of time and distance 
from the centre of impact by small steel cylinders and lead bullets. 

The time resolution of their measurements was only 50 ms, and they found 
that the temperature decayed rapidly as a function both of time and of 
distance away from the impact. Since in the current experiments the 
temperature at the centre of impact is observed within 1 ys of the shock 
wave arriving at the free surface, the high observed temperatures are not 
inconsistent with the maximum increase of 18°C observed 1.2 (or approxi- 
mately the crater radius) cm from the impact centre by Schneider and Stilp. 

One obvious omission from the Walsh and Christian formalism is the 
effect of elasto-plastic work. Although the metals have low Hugoniot 
elastic limits, they retain some rigidity after yielding, and may undergo 
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stress-hardening. The latter was reported by Fowles (1961), who demon- 
strated that an elasto-plastic equation of state should be used for 
aluminium shocked to pressures up to ^15 GPa. This concept was developed 
in detail by Lee and Liu (1967) and Lee and Wierzbicki (1967) who developed 
a treatment allowing direct calculation of temperature along the Hugoniot 
curve which incorporates the effects of material strength and finite 
uniso tropic strain. Foltz and Grace (1969) carried out the analysis for poly- 
crystalline aluminium having a yield strength comparable to that for 
aluminium- 20 24; whilst their analysis may not be strictly valid for the 
alloy used in these experiments, a comparison of their values for 
Hugoniot temperature (the dash-dot line in Figure 10(h)) and the measured 
residual temperatures is interesting. As can be seen, their calculated 
values are considerably in excess of the Walsh and Christian values at 
low pressures, but converge with them at higher pressures, which is 
precisely the behaviour observed in the measured residual temperatures. 
Although the release path is not certain, it has been proposed that for 
metals the release from shock pressure P occurs in two stages: first 
an elastic release (isentropic) to a pressure P„ - 2Y, where Y is the 

n 

Hugoniot elastic limit, followed by plastic release parallel to the 
Hugoniot (see e.g., Fowles, 1961; Al' tschuler, 1965). The post-shock 
temperatures indicated by the dotted curve in Figure 10(b) were derived 
from Foltz and Grace's Hugoniot temperatures assuming this form of 
release path with Y = 0.8 GPa. These values are in reasonable agreement 
with experimental observations (which may be slightly high due to changes 
in emlssivity). Unfortunately, a similar calculation cannot be made 
for an alloy like stainless steel where the constants required are not 
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really defined. However, the relationship between the observed and 
Walsh and Christian values (Figure 10(a)) is similar to that observed 
for aluminium, and plastic deformation of the lattice has been observed 
in recovery experiments for pressures up to 50 GPa (e.g. , Murr, 1975; 
Smith, 1958). The post-shock temperatures indicated by the solid line 
in Figure 10(a) were calculated from the Hugoniot temperatures assuming 
a simple isentropic release path; some allowance for elasto-plastic 
behaviour may be made by assuming a release path as indicated above and 
calculating the residual temperatures directly using 8. This yields 
the dotted line in Figure 10(a), hut the agreement between experiment 
and theory is not substantially improved. 

It is concluded that elasto-plastic work, which is not included 
in the Walsh and Christian formalism, causes significant heating at 
low pressures resulting in large differences between measured and 
calculated temperatures. However, at higher pressures (&30 GPa for 
aluminium., ^50 GPa for stainless steel) the Walsh and Christian approach 
appears to predict values close to those measured experimentally. 

c. Silicates 

A comparison between the observed and calculated post-shock temp- 
eratures for quartz is presented in Fig. 11. The agreement between the 
observations and the residual temperatures calculated directly from 
release adiabat data by Mashxmo et al. (1978) is very good, and this 
method of calculation is certainly to be preferred. However, the 
values calculated by Wackerle (1962) corrected for an initial temperature 
of 24°C are also remarkably close, to those observed except for the 
lowest and highest pressure points. The agreement may he improved by 
determining residual temperatures from the Hugoniot temperatures using 



22 


minimum release volumes estimated from the free surface and particle 
velocities given by Wackerle following the method described by Lyzenga 
and Ahrens (1978) or the release volumes given by Grady et al. (1974). 

Temperatures calculated using the method of Ahrens et al. (1969) 
using values of 6 and 8 GPa for the Hugoniot elastic limit (the lower 
and upper bounds given by Wackerle for z-cut quartz) were much too low; • 
the observed values are in fact greater than the calculated shock 
temperatures. Similar discrepancies were also noticeable for bronzite 
and forsterite, and may arise in the determination of the isotherm; 
it suggests this calculational technique is inadequate in the pressure 
range studied. 

For bronzite, none of the existing calculational tecnniques produced 
adequate agreement between theoretical and observed post-shock temperatures 
(Fig. 12) and since the observed values are generally in excess of the 
calculated shock temperatures some process must occur which results in 
greater heating on compression. (This may be related to the effects of 
sample porosity: the apparent temperature drop between 20.7 and 25.0 GPa for 
InSb measurements probably reflects the fact that the latter sample had a 
lower density of fine cracks than usual.) Attempts were made to improve the 
Walsh and Christian calculations by using a temperature, as well as volume, 
dependent Debye temperature in calculating the specific C via the Debye 
model. However, the post-shock temperatures are relatively insensitive to 
these refinements and no marked improvement in agreement was produced. The 
calculations of residual temperature are also fairly insensitive to the 
behaviour of Grtineisen's parameter provided it is the same on compression and 
on release. A model in which the Griineisen parameter is higher on compression 
than release yields higher residual temperatures. Although we cannot defend 
such a model in detail, tills Ldeo does in face have a good physical basis. 
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since the shock and 

release processes are certainly different. In particular, the observations 
of release paths lying beneath the Hugoniot (e.g., Pig. 1) and the 
apparent hysteresis in the shock and release process are not fully 
understood; they may' be related to the behaviour of y. A model in 
which Y 0 (=Y(V 0 ) was assumed to be 2.5 on compression and 1 on release 
produced somewhat better agreement, especially at high pressure* 
but a Grlineisens parameter of at least 3 during shock compression would 
be required to produce agreement at pressures less than 15 GPa. In this 
case the calculated release temperatures (assuming Yo = 1) would be too 
high at pressures of ^25 GPa, so a more complicated behaviour of y would 
have to be postulated. This does not seem justified since the theoretical 
Hugoniot does not adequately describe the effects of varying sample 
porosity, and the high temperatures may be largely due to this; further- 
more, the effect of the phase change has been neglected — although this 
would provide some justification for changing the behaviour of the 
effective y at 15 GPa. 

The results of the temperature calculations for forsterite are 
shown in Fig. 13; at this stage it should be noted that the 10 and 15 
GPa measurements using the InSb detector are probably contaminated by 
radiation from the sample driver Interface, as discussed previously. 

The observed temperatures are considerably in excess of the values cal- 
culated for the non-porous polycrystalline forsterite using the equili- 
brium Hugoniot (curves A, A T , Fig. 13(a)), and are also greater than the 
values calculated for the theoretical Hugoniot (C, C') although there is 
some indication that the measurements tend towards the latter at high 
pressures. The values calculated for the "porous" equilibrium Hugoniot 
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(B, B* ) are much higher than the observations except for the doubtful 
InSb points; this is not surprising since the samples were only £1% 
porous, and not 4%. As with br onsite, a better fit to the observations 
may be obtained using the Walsh and Christian approach but assuming that 
Yo is higher on compression than release; a fairly good fit resulted 
from using the theoretical Hugoniot with Yo ~ 2.5 on compression and 
Y 0 = 1 on release. However, although temperatures calculated by the 
conventional methods of Walsh and Christian and Ahrens et al. are much 
lower than the observations at low pressure, they converge at higher 
pressure, which is reminiscent of the behaviour observed in stainless 
steel and aluminium and ascribed to elasto-plastlc effects. A detailed 
investigation of the process of dynamic yielding in forsterite would 
indicate whether such effects were possible here; if the behaviour of 
forsterite is indeed similar to that of metals, then the release path 
may be approximated by the two stage elasto-plastic path that was 
described earlier and is sketched in fig. 14. In this case an estimate 
of the release temperature may be obtained by calculating the energy 
difference between the shock and release paths (the shaded area in 
Fig. 14) and equating it to the increase in internal energy as in 
equation 8. This results ia. the temperatures shown by the heavy line 
in Fig. 13(b); the agreement between these values and the observations 
is good at ^10 GPa, but large differences exist at higher pressures, as 
would be expected if the material underwent a gradual loss of strength 
that was complete by ^25 GPa. (Experimental evidence indicates that the 
loss of strength is in fact complete by 'uZO GPa; see Appendix B.) This 
calculation ignored the Fraction oF energy going into permanent deformation 
of the lattice, which may be estimated as follows: the free energy of 
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3 

a crystal is increased by about pb /2 per atomic length of dislocation, 
where p is the shear modulus and b is the average length of the Burgers 

O 

vectors of the crystals. With b^5A, p - 80 GPa, and a dislocation 
Xl -2 

density of 10 cm ■ , which is probably rather high for these pressures 
(R. Jeanloz, personal communication, 1978), allowing for the deformational 
energy results in a 15% or smaller reduction in the increase in thermal 
energy due to shock compression and release. This is equivalent to a 

(or smaller) decrease in release temperatures from the values plotted 
in Fig. 13(b). 

Discussion 

The experimental technique described here has proved useful for 
measuring post-shock (brightness) temperatures in a variety of metals 
and silicates. Initial experiments produced results which are, in general, 
consistent at two wavelength ranges (from 4.5 to 5.75p and from 7 to 
14p) and for different samples of the same material shocked to similar 
pressures. Although the consistency of the results suggests that the 
measured temperatures are Indeed representative of the residual temperatures 
in the shocked samples, uncertainties still exist in the precision of 
the metals because of the lack of information on the effect of surface 
process, which are likely to he important in metals, and changes In 
emissivity under shock conditions. In particular, little data is 
available on the emissivity of the silicates studied, even under normal 
conditions; and if such data were available more accurate correction 
factors could be estimated, although the possibilities of non-equilibrium 
radiation with an emissivity greater than unity and changes in emissivity 
under shock conditions would remain. Experimental investigation of the 



26 


latter, if possible, would be most valuable. One effect that may 
influence the measurements of residual temperatures' in aluminium is the 
ejection of material from the surface; this has been studied by Asay 
et al. (1976) for pressures of ^25 GPa and is quite significant. 

A further study by Asay (1977) showed that material ejection was highly 
dependent on the rise-time of the shock wave, so the effect of material 
ejection on temperature measurement might be investigated by determining 
residual temperatures for different shock rise- times. 

Measured temperatures in quartz, which is known to undergo a 
complete loss of strength after the Hugoniot elastic limit is exceeded, 
were in good agreement with those calculated by Mashimo et al. (1978) 
from release adiabat data, and in reasonable agreement with those derived 
from the equilibrium Hugoniot of Wackerle (1962) using estimated release 
volumes. These calculations did not include the affects of the quartz 
to stishovite phase change which begins at VL5 GPa; however there was 
only a slight change in the rate of increase of measured temperature with 
pressure in this region. 

For both the metals studied, and for forsterite and bronzite, the 
measured temperatures were considerably in excess of those predicted using 
conventional fluid rheology assumptions, although the absence of release 
adiabat data for forsterite and bronzite meant tnat a comparison with 
temperatures calculated in the optimal manner (i.e. using equation 8) 
was not possible. In all these cases some process is required which 
leads to higher calculated shock temperatures, especially at low pressures, 
and less energy loss on release. For metals, this is probably elasta- 
plastic work, which dominates at low pressures: reasonable agreement was 
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found between the elasto-plastic calculations of Foltz and Grace (1969) 
and measurements for aluminium. An attempt was made to model such a 
process in bronzite and forsterite by allowing Griineisen's parameter to 
be higher on compression than on release, and for forsterite reasonable 
agreement was obtained with a 2.5 on compression, y Q = 1 on release 
(note that tinder standard conditions the thermodynamic Griineisen' s parameter 
is 1.17). No strict physical justification for postulating a different y 
for shock compression from that on unloading is attempted here; however, 
this type of behaviour does provide better agreement between theory and 
observations although the model is somewhat non-unique. The differences 
between observed and predicted release adiabats and the hysteresis observed 
in shock unloading suggest that some change in material properties may occur 
prior to release from the shock state, and thus irreversible behaviour of the 
enharmonic properties is not incredible. . Another possibility is that elasto- 
plastic work may also be important: although quartz has been observed to 
undergo a complete loss of material strength beyond the Hugoniot elastic 
limit, the same has not been demonstrated for forsterite. Calculations of 
residual temperatures in which it was assumed that forsterite behaved in 
an elasto-plastic manner similar to that proposed for metals yielded values 
that agreed well with observations at 10 GFa but were much higher at 
higher pressures. This suggests that forsterite undergoes a much more 
gradual loss of material strength beyond its elastic limit than does 
quartz, and that elasto-plastic effects may be important at pressures 
below *u25 GPa. 

Classifications of shock metamorphism such as that of Staffer, 

(1971) are based on a detailed analysis of the physical characteristics 
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of the specimen (such as the amount of fracturing, undulatory extinction 
or presence of glass.) Estimates of the temperature range corresponding 
to each class or metamorphism have then been made using calculated values 
of the shock and post-shock temperatures: Stoffler used the values of 
Wackerle (1962) for quartz and Ahrens et al. (1969) for feldspars. The 
results of the present study indicate that the former are fairly reliable 
below ^20 GPa (when corrected for the right initial temperature) , hut 
the ealculational technique of Ahrens et al. yielded temperatures that 
were substantially lower than those observed, so that the values used 
by Stuffier for feldspars are also likely to be low. Indeed, observations 
of recovered shocked samples suggest that this is the case both in 
feldspars (e.g. Gibbons, 1974) and dunite (e.g. Reiraold and Stdffler, 

1978). Measurements of the electrical conductivity of shocked dunites 
(Schulein et al. 1978) show increases that are also more easily explained 
if the actual temperatures are higher than those calculated. The main 
implication of this study for impact metamorphism is that temperatures in 
forsterite and bronzite appear considerably higher than the theoretical 
values, especially for porous samples, and so the effects of impact 
metamorphism on basic rocks such as lunar basalt may differ from those 
expected on the basis of theory and lead to erroneous conclusions 
regarding the nature of the impact. Since the partitioning of energy 
during impact processes is important in calculating the thermal history 
of an accreting planet, residual temperature data such as those pre- 
sented here could provide useful constraints on the temperature distribution 
during accretion. 


One important assumption used in the temperature calculations was 
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that equilibrium existed behind the shock-front. Temperatures measured 
in these experiments were generally higher than those calculated, but 
other factors leading to errors in the calculation make it hard to assess 
the validity of this assumption. A systematic investigation of differences 
between shock effects seen in shock recovery experiments and in naturally 
shocked samples where the peak pressures are of longer duration could 
help clarify this point. Such differences certainly exist - for example, 
the substantial amount of coesite observed in natural, but not shock- 
recovery, specimens (StSffler, 1972) and might yield important information 
on the nature of equilibrium in shock processes. 

It has been emphasised that the release path is of critical importance 
to the dete rmina te on of post-shock temperatures, yet it has generally been 
assumed that the release Is isen tropic; if this is not the case, then 
discrepancies between measurement and theory are not surprising. Recently 
Kieffer and Delany (1979) have placed new and useful constraints on the 
assumptions of isentropic flow in solids under decompression. Viscous 
heating on release would lead to higher residual temperatures. 


9 

Release adiabats which lie below the Hugoniot (see Fig. 1) lead 
to higher residual temperatures in better agreement with the observations, 
yet it seems strange that if the material is essentially unchanged by 
shock compression and release, it should occupy a smaller volume on 
release when it is hotter than it was initially. Measurements of the 
parameters on release adiabats are often made at relatively high pressures 
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(P >P_ »0, where P„ is the shock pressure and P„ the point at which 
release measurements are made) and so yield in£ormation on the initial 
slope of the release adiabat which may then be extrapolated to yield the 
release volume. Such an extrapolation is illustrated schematically in 
Pig. 15 by the dashed line: this yields a release volume smaller 
than Vq 5 however the actual path may be as shown by the dotted line. 

Since the area between the Hugoniot and the two release curves is the 
same, they will give the same residual temperatures if (8) is used. 

Assuming that the material is unchanged, then the release volume may 
be estimated from the post-shock temperature and the mean thermal expansion 
coefficient a^, and can provide additional constraints on the release 
path. Table 10 lists the estimated release volumes for the materials 
studied: in fact, for the silicates and stainless steel they are not 
very different from the initial volumes. 

It was hoped that this study would yield some definite information 
on the behaviour of Grimeisen’ s parameter, which is critical for the 
reduction of shock wave data to the form needed for comparison with 
density-depth profiles within the earth. However, post shock temperatures 
are not very sensitive to the behaviour of y (unless it is different on 
shock compression and release) and appear more greatly influenced by the 
release path, in particular the release volume. There seems to be no 
need to postulate unusual behaviour of y for quartz to V25 GPa, although 
the high observed values of residual temperatures in forsterite and 
bronsite suggest that present rheological models used in shock temperature 
calculations, at least below '^25 GPa, are either inadequate, or the 
effective Griineisen parameter is high (y Q ^2. 5) on compression and 
decreases to near its zero pressure value (y^l) upon adiabatic release. 
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Table 1 

CALCULATED POST SHOCK TEMPERATURES 
IN FUSED QUARTZ AND OLIGOCLASE, 

using the Mie-Gruneisen Theory (A) or Release Adiabat Data (B) 


Post Shock Temperatures, °C 


Shock Pressure 

Fused Quartz 

Olxgoclase 

GPa 

A 1 

B 2 

A 3 

B 4 

10.0 

0 

80 



15.0 

0 

450 



18.0 



27-35 

269-386 

25.0 

0 

1220 



27.2 



129-206 

> 742 

30.0 

470 

1480 



40.0 

1860 

2180 



41.7 



327-395 

>1031 

50.0 

3310 

2820 




1. Wackerle (1962) 

2. Gibbons and Ahrens (1971) 

3. Ahrens et al . (1969) 

4-. Ahrens and O’Keefe (1972) 


Table 2 


ESTIMATED SIGNAL CORRECTION FACTORS 

Correction factor — y (values to nearest 0.05) 

Black Body signal 



T = 

400°K 

T “ 

600°K 


InSb 

HgCdTe 

InSb 

HgCdTe 

Quartz (polished) 

0.90 

0.80 

0.85 

0.80 

Quartz (rough) 


0.90 


0.90 

sio 2 

0.90 


0.90 


Dunite (polished) 

0.85 

0.90 

0.80 

0.90 

(Bronzite 

0.9 

0.85 

0.9 

0.85) 


Table 3 


POWER LAW PITS TO CALIBRATION CURVES 

t) 2 , 

S = a(T~24) D x 10 ; coefficient of determination r 

Black Body Stainless Steel Alumininm-2024 



a 

b 

2 

r 

a 

b 

2 

r 

a 

b 

2 

r 

InSb 

7.21 

1.87 

.99 

2.21 

1.93 

.98 

3.08 

1.71 

.95 

HgCdTe 

10.9 

1.45 

.97 

10.05 

1.32 

.97 

2.68 

1.43 

.98 
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Table 4 

MEASURED POST-SHOCK TEMPERATURES 
IN STAINLESS STEEL-304 
Temperature, °C 



InSb 


Hf 

*CdTe 

Peak 

| 

Pressure, GPa 

SS 

BB 

SS 

BB 

(InSb) 

. ^ 

( 

i 

11.5 



80 

60 


ik 

4 

11.7 

no 

75 



250 


13.0 

125 

80 



600 

! 

14.5 



130 

85 



16.0 

145 

100 

145 

95 

830 


23.0 

195 

130 



1530 


24.2 



200 

130 


T 

] 

43.0 

325 

230 




1 

1 

i 

50.0 

355 

250 



1820 

| 


SS - calibration using stainless steel 
BE = black body temperature 

Uncertainties in temperature: ±15° below 150° , ±10° above 150°. 
T = 24° C 
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Table 5 

MEASURED POST-SHOCK TEMPERATURES 
IN ALUMINIUM-2024 
Temperature, °C 



InSb 


HgCdTe 


Peak 

jssure, GPa 

AL 

BB 

AL 

BB 

(InSb) 

10.5 

125 

50 




11.5 



140 

65 


12.5 

135 

55 



1250 

15.0 

150 

60 




15.7 



155 

70 


18.5 

175 

75 

185 

80 

1430 

25.0 



220 

90 


27.0 

230 

105 



2200 

32.5 

250 

120 




33.0 

260 

127 



3800 


AL = calibration with aluminium 

BB = equivalent black body temperature 

Estimated uncertainties: ±20°C below 200°C, 


T = 24°C 


±10°C above 200°C 
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Table 6 

MEASURED POST-SHOCK TEMPERATURES IN QUARTZ 
Temperature, °C 


IuSb HgCdTe 


Pressure, 

GPa Flash 

BB 

Corr. 

Flash 

BB 

Corr, 

5.0 

235 

80 

87 




5.5 




110 

. 75 

85 

8.0 

180 

100 

105 




9.5 

225 

120 

125 

320 

115 

127 

10.8 

245 

155 

162 




11.5 

252 

160 

170 




15.0 

340 

177 

187 




15.5* 




— 

160 

180 

17.5 

377 

185 

195 




19.5 




706 

320 

340 

20.0 

390 

242 

255 




21.5 

425 

250 

265 




*This was 

a very faint 

record. 

and may 

not be reliable 

■ 


Uncertainties: ±10°C below 100°C,±5°C above 100°C 


T - 24°C 
o 
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Table 7 

MEASURED POST-SHOCK TEMPERATURES IN FORSTERITE 
Temperature, '“C 


InSb HgCdTe 


Pressure, GPa 

Flash 

BB 

Corr. 

BB 

Corr. 

7.5 

180 

<50° 

(no detectable 

rise) 


9.6* 

237 

105 

115 

65 

77 

15.0** 

260 

136 

145 



18.0 




105 

112 

20.2 

285 

140 

152 



21.0 




120 

125 

24.0 




160 

165 

24.5 

270 

148 

160 



28.0 

300 

156 

167 




*Twa peaks (175, 237° C); residual temperature corresponds to the 
difference between the levels after the second and first peaks. 


**Two peaks (135, 260°C); residual temperature estimated as before. 


Estimated uncertainties: ±10° C below 100°C, ±5°C above 100° C 


T = 24° C 
o 
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Table 8 

MEASURED POST-SHOCK TEMPERATURES IK BAMBLE BR.ONZITE 
Temperature , °C 


Pressure, GPa Peak BB Corr. BB Corr . 


6,0 


< 50° 

(no detectable rise) 


10.3 

123 

100 

105 


11.0 



no 

120 

14.8 



145 

160 

15.5 

185 

147 

157 


20.7 

225 

200 

213 


21.5 



185 

200 

25.0 

200 

175 

185 


26.0 



225 

240 

Uncertainties : 

±10 D C below 

100°C, 

±5°C above 100°C 



T = 24°C 
o 
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Table 9 

SOME CONSTANTS RELEVANT TO THE CALCULATION 
OF TEMPERATURES IN SHOCKED SILICATES 



p” 1 

cm^/gm 

Y P 

Od 1 

°K 

m 

K 

os 

GPa 

K* 

os 

Quartz 

.377 

. 703 2 

1050 

20.03 

2 

37.7 

6.4 2 

Forsterite 

.310 

2 3 
1.17 Z,J 

900 

20.12 

(126.7“*, 

\128.8 4 

5.37 3 

Bronzite 

.298 5 

.907 

950 

20.96 

103. 5 5 

9.59 5 


.307 

5 

1.56 3 



105. 0 7 

5.3 6 


Y 0 = thermodynamic Gruneisen parameter 
0 q - Debye temperature 

m = mean atomic weight 

K = zero pressure adiabatic bulk modulus 
os 



1. Debye temperatures derived from fitting specific heat data from 
J.A.N.A.F. Tables. 

2. Values from Anderson et_ al . (1968) 

3. Kumazawa and Anderson (1969) 

4. Graham and Barsch (1969) 

5. Frisillo and Barsch (1971) [Mg^ gFe^ 2^ Si ^*3^ 

6. Chung (1971); preferred value yielding better fit to Hugoniot data. 

7. Kumazawa (1969) 

Note: 7 

v was generally assumed constant. 


Table 10 


ESTIMATED RELEASE VOLUMES 



Al-2024 

SS-304 

Quartz 

Forsterite 

Bronzite 

v 0 (cm 3 /g) 

0.359 

0.127 

0.377 

0.310 

G.302 + 

% 

8 x 10“ 5 

5 x 10“ 5 

3.5 x 10* 

-5 -5 

3 2.6 xlO 3 

2.4 x 10" 

Pressure, GPa 






5 

— 

— 

0.3776 

^0.310^ 

*0.3 02** 

10 

0,361 

0.1274 

0.3782 

O'. 3104 

0.3026 

15 

0.363 

0.1276 

0.3789 

0.3105 

0.3029 

20 

0.364 

0.1279 

0.3799 

0.3107 

0.3032 

25 

0.365 

0.1281 

— 

0 . 3111 

0.3031/ 

0.3035 


+ Theoretical zero-pressure density, (Ahrens & Gaffney, 1971) 
Actual initial volume *0.304 cm ‘Vg 

-H- Mo detectable temperature change 
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Appendix A Detector parameters. 

Both detectors were purchased from the Santa Barbara Research Centre, 
Goleta, California. 

a) InSb 

The InSb detector used was a circular chip 1 mm. in diameter 
having a detectivity of 5 x 10^ cm Hz^^/watt when operated at 77°K. 

It was operated with a fast matched pre-amplifier (S.B.R.C. Model A-230) 
consisting of a current mode operational amplifier with a feedback 
resistance of 1 kft and a non-inverting voltage mode post-amplifier; 
this stage had a gain of 500, and upper and lower 3db frequencies of 
20 Mi-Iz and 1.35 kHz respectively. For use in measuring post-shock 
temperatures an additional amplifier with variable gain (from 1000 to 
30,000) was used (Raikes, 1978). The minimum system rise time of 
^O.lps is controlled by the detector chip itself. Sapphire windows were 
used with this detector. 

b) HgCdTe 

-2 2 

The chip used had an area of ° x 10 cm , a detectivity of 
9 1/2 

6.94 x 10 cm Hz /watt and a rise time of 100 nsec. It was used with 
a matched amplifier, S.li.R.C- Model A-120, having a gain of 1000. The 
amplifier consisted of an a-c coupled voltage mode amplifier plus a 499 
load resistor and circuitry to produce the bias current of 10 mA required 
by the detector; its upper and lower 3db frequencies were 10 MHz and 50 Hz 
respectively. The rise time of the detector-amplifier system is ^ 0.05 ps 
however, for operation at low signal levels it was found to produce an 
unacceptable level of very high frequency noise, and so had to be operated 
with a filter which raised the rise-LLme to 'v 0.75ns. Barium fluoride 
or Irtran-2 (Kodak) windows wore used with the HgCdTe detector. 
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Appendix B Dynamic Yielding in Forsterite 

A detailed study of the dynamic yielding of single crystal forsterite 
has not been carried out. For the purpose of the present work a single 
measurement of the Hugoniot elastic limit along the c-axis of forsterite 
was undertaken in order to place some bounds on the rheological behaviour 
at relatively low shock stress. 

Shock loading and recording were carried out with a modification of 
the 40 mm gun system described in .Ahrens at al. , 1971, 1973, although the 
optical recording was carried out using a Model 339B Beckman and Whitley 
continuous writing streak camera with a recently constructed xenon light 
source (described in Goto et_ al . , 1978) . 

The sample, machined from the same aliquot as that utilised for 
the post-shock temperature measurements, was mounted on a 1.5 mm thick 
aluminium-2024 driver plate and impacted with a 4mm thick aluminium-2024 
flyer plate. The pertinent measurements obtained in a single experiment 
(Figure B-l, Table B-l) provide a rather unequivocal measurement of the 
Hugoniot elastic limit amplitude, 8.68 ± 0.67 GPa, for a 3.7 mm thick sample. 
The velocity of the first shock, 8.72 km/s, is close to, but slightly higher 
than, the values of 8.543 km/s (Kumazawa and Anderson, 1969) or 8.564 km/s 
(Graham and Barsch, 1969) measured for the longitudinal velocity in the 
(3) direction. A similar relationship between the two velocities has been 
reported for quartz in various orientations (Wackerle, 1962; Fowles, 1967). 
The theoretical relation of the higher order elastic constants to the 
velocity of finite strength elastic waves, although reported for hexagonal 
symmetry (Fowles, 1967) has not been worked out for orthorhombic crystals. 
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The arrival time of the second shock, as indicated by the marked change 
in slope observed in the image of the inclined mirror (Figure B-l) , can 
be interpreted in several ways depending on the rheology assumed for the 
sample material which has been compressed, and subsequently released, by 
the first shock (Wackerle, 1962) . Our preferred analysis is simply to 
assume that the second shock velocity is given by 



d + U l fs ^2 V 


B-l 


where d is the sample thickness, U^ s the free surface velocity and 
(t^-t^) and (t^-tQ) are the time interval between the arrival of the 
first and final shock and the travel time of the second shock to the free 
surface respectively (Ahrens et_ al , , 1968) . Equation B-l ignores the 
interaction of the reflection of the first shock with the oncoming final 
shock front. The parameters of the fLnal shock state (Table B-l) were 
obtained from an impedance match solution using the equation of state 
parameters given by McQueen et al. (1970) for aluminium-2024. Equation 
B-l yields a final shock state which agrees well with the pressure-density 
trajectory calculated via the Murnaghan equation using parameters obtained 
from the ultrasonic data. If, on the other hand, the expression for the 
final shock velocity given by Ahrens et al . (1973) , which assumes an elastic 
interaction of the reflected elastic shock with the oncoming final shock, 
is used a state with a greater compression is calculated (Figure B-2) . 

Although the difference between Lhe two solutions is not large, the calculation 
assuming no interaction, implying a loss of strength at pressures corresponding 
to the final shock state, appears to be definitely more consistent with 
the ultrasonic data. The relatively high Hugoniot elastic limit, followed 
by a (gradual) loss of shear strength, is consistent with the rheological 
behaviour inferred from the post-shock temperature measurements (Table 7 ). 


Table B-l 


Equation of State Data for Forsterite Shocked 
Along the (001) Direction 
(Shot 431) 


Sample Mass 
Bulk Density 
Projectile Velocity 


3.1546 ±0.0001g 
3.2158 ±0.005 gm/cm 3 
1.668 ±0.002 km/s 


Shock velocity 
Shock density 
Shock pressure 

Free surface 
velocity 


First Shock State 

8,717 ±0.068 km/s 

3.334 ±0.008 gm/cm 3 

8.68 ±0.67 GPa 

(Hugoniot elastic limit) 

0.619 ±0.042 km/s 


Final Shock State 
6.554 ±0.067 km/s 
3.5716 ±0.0041 gm/cm 
17.32 ±0.05 GPa 


Particle velocity 


0.7246 ±0.0022 km/s 
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Figure Captions 

Figure 1. Schematic diagram of Hugoniot and possible release adiabats. 

The solid release curves lie above the Hugoniot, and the 
dashed ones below it. 

Figure 2. Schematic diagram of the experimental configuration. An 

oscilloscope recording the detector output at 50 ps/ div is 
triggered by the passage of the projectile past the first 
laser beam. A record having greater time resolution (5 ps/div) 


Figure 3. 

j 

i 


Figure 4. 


Figure 5. 


i 

i 

is obtained from an oscilloscope triggered by the impact of the \ 

| 

projectile with the shorting target. 

(a) Response curve for InSb detector material. 

> 

(b) Response curve for HgCdTe detector material. 

D*(A,v) is Lite area independent sensitivity at frequency v 
and wavelength A. 

Shaded areas indicate typical variations in detector sensitivity. j 

(Details from the Santa Barbara Research Centre Catalogue) < 

Infra-red transmission scans for the silicate materials studied. j 

■i 

The operating range of the filtered InSb detector and the start 

£ 

of the HgCdTe band are also shown. l 

Calibration curves for the InSb detector, operating in the $ 

wavelength range 4.5 to 5.75p. The open and solid circles in 

'1 

the aluminium curve are the results of two different calibration runs. 


Figure 6. 
Figure 7- 

Figure 8. 


Calibration curves for the HgCdTe detector. 

InSb detector output records for aluminium and stainless steel. 
Tj = residual temperature, = peak temperature, A = air shock. 
InSb detector output records for the silicates studied. 

T| = flash temperature, (or T) = residual temperature. 


A = air shock. 
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Figure 9. HgCdTe detector output records for the silicates studied. 

T — residual temperature. 

Figure 10. Comparison of observed and calculated temperatures in 

(a) stai nl ess steel and (b) aluminium. The dashed curves are 
the Hugoniot temperatures calculated using the Walsh and Christian 
technique, and the solid curves the corresponding release temperatures. 
The dotted curve in (a) represents residual temperatures calculated 
directly as su min g elasto— plastic release. The dot-dash and dotted 
curves in (b) are respectively the shock and post-shock temperatures 
obtained using the treatment of Foltz and Grace (1969) . 

Figure 11. Comparison of observed and calculated temperatures in quartz. 

In (a) the dashed line and solid line are the values of Wackerle 
(1962), corrected for T ~ 24°C. Asterisks (*) indicate residual 
temperatures determined from the plotted Hugoniot temperatures 
using estimated release volumes. In (b) the dashed lines are 
Hugoniot temperatures calculated as per Ahrens et al. (1969) for 
Hugo ni ot elastic limits as indicated; the solid line is the release 
temperature for an H.E.L. of 6 GPa. The heavy solid line shows 
the values of Mashimo et al. (1978). 

Figure 12. Observed and calculated temperatures for Bamble bronzite. Solid 
curves are post-shock temperatures, broken curves Hugoniot 
temperatures. 

(a) A, A': theoretical Hugoniot y =* .907 

B, B 1 : theoretical Hugoniot, y = 2.5 on compression, 

1 on release 

C, C T : actual Hugoniot, y - .907, Hugoniot elastic limit 

of 6.7 GPa (Ahrens and Gaffney, 1971). 
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(b) D, D' : theoretical Hugoniot, y = 1.57 

E, E’: actual Hugoniot, y ~ 1.57, Hugoniot elastic 
limit of 6.7 GPa 

Figure 13. Observed and theoretic il temperatures for forsterite 

(a) A, A 1 : polycrystalline forsterite, equilibrium Hugoniot 
B, B’ : ^4% porous forsterite, equilibrium Hugoniot 

C p C 1 : theoretical Hugoniot 

(b) D, D ? : theoretical Hugoniot, y - 2.5 on compression, 

y = 1 on release 

E, E’: actual Hugoniot, assumed Hugoniot elastic limit 
5 GPa 


F : actual Hugoniot elastic limit 9 GPa shock temperatures 

Note: The measured H.E.L. was 8.7 GPa. 

Figure 14. Schematic diagram illustrating elasto-plastic release path. 

Figure 15. Two release paths having the same initial slope and area 
underneath yet different final volumes. 

Figure B— 1 Static and dynamic streak images produced upon shock compression 
of c-cut single crystal forsterite to 17 GPa. (a) Static image 
as seen through streak camera, (b) Streak image demonstrating 
two-wave shock structure recorded by inclined mirror. The angle 
Y, is related to the value of the Hugoniot elastic limit, see 
e.g. Ahrens et al. (1968). 


Figure B-2 Shock pressure versu . density plot for forsterite showing the 
position of the Hugoniot elastic limit. The isentropic 
compression curve was calculated using a value of the bulk 
modulus and its pressure derivative of 128.8 GPa and 5.37 
(see Table 9). 
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